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INTRODUCTION  -  Tasks  Undertaken 

In  the  first  and  second  years  of  contract  support,  we  have  followed  the  assi^ed  tasks  on 
the  imaging  capab^ty  and  tumor  characterization  in  human  breaste  using  TRS  (time  resolved 
spectroscopy)  for  1500  measurements.  We  have  computed  optical  characterization;  i.e.,  Pa 
(absorption)  and  m’  (scattering)  measiuements  for  20,  34,  and  6  patients  in  the  first,  second  and 
third  years,  respectively.  The  statistics  of  the  60  patients  with  various  stages  of  tumors  ate  available 
and  compared  to  17  healthy  normal  controls.  For  the  imaging  capability,  we  have  shown  that  with 
the  current  technologies  and  the  current  limitations,  we  can  detect  objects  of  ~  1  cm.  Longer  times 
and  transmission  studies  will  improve  the  resolution. 

In  the  past  year  (year  3)  we  have  studied  intensively  novel  modalities.  The  specified  tasks 
required  360  studies.  1)  TRS  or  PMS  0>hase  modulation  spectroscopy)  can  be  used  to  iden^ 
oxygenation  and  scattering  coefficient  of  tumors.  2)  TRS  may  not  be  the  best  method  to  provide 
imaging  mainly  because  of  time  requirements.  Thus,  we  concluded  that  the  TRS  is^  perhaps  more 
suitable  for  accurate  measurements  of  the  optical  characterization,  such  as  hwnoglobin  oxygenation 
in  the  tumor  and  scattering  coefficients.  In  addition,  our  interests  in  the  optical  imaging  now  include 
the  application  of  PMS  and  CWS  (continuous  wave  spectroscopy). 

In  the  following  sections,  we  report  on  various  tasks  as  follows: 

1)  Optical  breast  imaging  using  TRS; 

2)  TRS  Study  of  qitical  characterization; 

3)  Present  and  future  plans  for  the  TRS  system  for  optical  charactoization  in 

tumor  bearing  breasts; 

4)  CWS  breast  imaging  with  contrast  agent; 

5)  Imaging  using  the  frequency  domain  method. 

1.  OPTICAL  BREAST  IMAGING  USING  TRS  (Years  1  and  2) 


1.1.  Introduction 

Despite  the  long  history  of  the  optical  imaging  technology,  it  has  not  been  commonly  used 
for  breast  cancer  detection.  Photon  diffusion  into  the  tissue  m^es  a  shadow  of  the  object  too 
obscure  to  use  for  imaging.  The  shadow  imaging  with  continuous  light  (diaphanography)  has  been 
used  in  Asia,  where  average  breast  size  is  relatively  smaller,  but  the  shadow  techiuque  does  not 
allow  real  imaging.  Today,  with  new  technologies,  such  as  the  time  gated  technique  (1),  more 
diffusive  photons  are  eliminated  and  therefore  better  images  can  be  obtained.  We  have  studied  the 
ability  of  the  time  gated  optical  technique  to  human  breast  cancer. 

We  have  chosen  to  develop  the  time  gated  optical  technique  for  breast  cancer  detectiwi  for  the 
following  reasons.  First,  the  breast  tissue  consists  of  a  relatively  low  scattering  media,  namely  fat, 
whose  scatter coefficient  is  as  low  as  3  cm’*.  Second,  the  breast  is  located  on  the  surface  of  the 
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body  and  is  accessible  to  ihe  optical  instruments.  Thirdly,  the  time  gated  technique  will  allow 
selection  of  photons  which  have  not  diffused  extensively  (die  extensively  diffused  photons  are  the 
ones  which  create  obscure  images).  We  selected  15%  of  total  photons,  which  diffused  relatively  less 
and  arrived  earlier  than  85%  of  rest  of  photons.  We  emphasize  not  only  possible  improvement  of 
sensitivity  of  detection  by  this  time  gated  technique,  but  also  high  specificity  of  cancer  detection  by 
functional  imaging.  It  creates  an  oxygen  image,  a  blood  volume  image  and  a  scattering  image,  and 
therefore  can  te  helpfiil  in  interpreting  the  cell  function  of  lesions.  Preliminary  studies  using  TRS  for 
characterizing  tumors  in  the  breast  have  been  published  elsewhere  (2). 

1.2  Subjects  and  Methods 

Thirty  five  women  who  had  breast  tumors  detected  by  X-ray  mammogram  were  subjects  for 
the  optical  imaging.  Hie  age  of  the  women  ranged  from  30  to  84  years  old.  Of  them,  15  patients  had 
had  ^  lesions  removed  (lumpectomy  operation).  20  patients  h^  a  single  suspicious  small  tumor 
shadow  in  the  mammography.  The  tumor  sizes  ranged  from  3  mm  to  1  cm  diameter.  The  post¬ 
operative  scar  tissue  contained  seroma.  The  size  of  the  post-operative  lesions  were  2  cm  to  5  cm  in 
diameter.  Approximately  thirty  percent  of  the  patients  h^  fibrocystic  disease. 

For  the  optical  measurement,  the  breasts  were  con^ressed  gently  with  two  surface  holders 
located  parallel  to  each  other  (Figure  1).  One  surface  holds  an  incident  light^ide  and  the  other  holds 
a  lightguide  coupled  to  a  detector,  photomultiplier.  In  between  the  two  lightguides,  time  resolved 
spectra  of  the  breast  are  acquired  through  transmittance  mode.  Then,  points  of  tire  lightguides  on  the 
two  surfaces  are  n»ved  5  rrrm  equally  to  obtain  a  two  dimensional  transmittance  image.  During  the 
acquisition,  lightguide  separation  is  fixed  and  constant. 


CHEST 


Fig.  1.  The  scheme  shows  the  set  up  for 
the  optical  breast  imaging.  The  breast  is 
gently  confessed  with  two  surface 
holders  located  parallel  to  each  other.  One 
surface  holder  contains  an  incident  light¬ 
guide  and  the  other  holder  contains  a 
lightguide  coupled  to  a  detector.  In 
between  the  two  lightguides,  TRS  data  of 
the  breast  are  acquit  through  trans¬ 
mittance  mode. 


Transmittance 


The  time  resolved  spectrometer  (TRS)  consists  of  two  wavelength  pulsed  lasers,  which  emits 
light  at  780  and  830  nm,  a  photomultiplier  (gallium  arsenide),  a  time  amplitude  converter  (TAG),  a 
photon  counting  system,  and  the  data  stored  in  a  personal  computer  through  an  A/D  converter,  lire 
time  width  of  Ae  puls^  laser  light  is  300  picoseconds  through  Ae  TRS  system  and  is  used  to 
detemmie  mstrument  fiinction.  This  time  resolution  of  Ae  mstrument  function  is  good  enough  for 
human  breast  spectra  to  be  used  wiAout  deconvolution  because  Ae  tissue  spectra  gained  Ae  half 
widA  more  Aan  5  times  larger  Aan  Ae  instrament  function.  The  two  tissue  spectra  at  Ae  two 
wavelengAs  are  acquired  simultaneously  by  time  sharing  (delaying  one  pulse  by  10  ns). 
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The  two  tissue  spectra  with  wavelengths  at  780  ran  and  830  ran  are  analyzed  to  obtain 
absorption  and  scattering  coefficients  with  a  curve  fitting  program,  which  in  principle,  is  based  u^n 
a  diffusion  equation  using  semi-infinite  boundary  conditions  (Patterson,  et  al)  (3).  Here,  scattering 
coefficients  are  mainly  calculated  at  maximum  photon  arrival  time.  This  requires  15  to  20%  of  eariy 
arrival  photons.  The  absorption  coefficient  is  mainly  calculated  by  the  decay  curve  of  the  spectra, 
requiring  25  to  50%  of  the  early  arrival  photons.  Hemoglobin  saturation  with  oxygen  (Sa02)  is 
calculatai  using  the  two  absorption  coefficients  (4).  As  a  presentation  of  blood  volume,  the  two 
absorption  coefficiraits  arc  averaged.  Scattering  coefficients  at  two  wavelengths  are  also  averaged  to 
yield  a  scattering  image. 

Three  independent  pieces  of  information  arc  used  to  process  three  images  from  a  set  of 
measurements:  1)  the  oxygen  image  is  based  upn  hemoglobin  saturation  with  oxygen,  2)  blood 
volume  image  is  from  absorptions,  and  3)  scattering  image  is  from  scattering  coefficients.  Chice  these 
three  are  calculated,  diey  are  fed  into  the  image  processing  program,  which  translates  the  location  of 
lightguides  on  the  breast  surface  as  a  pkel  of  the  image.  It  therefore  creates  images  of  a  5  mm^ 
resolution.  Usually  a  5  x  7  pixel  image  is  acquired. 

1.3  Results 


It  is  found  that  the  ability  to  detect  breast  lesions  by  the  time  gated  optical  method  depends 
upon  the  size  of  lesions,  the  thickness  of  the  breast  and  contents  of  breast  tissue  (fat  and  fibrous 
tissue).  Usually  post-operative  lesions  can  be  seen  in  our  optical  images,  but  a  small  cancer  mass  of 
less  dian  1  cm  in  dameter  is  not  detectable  with  5  to  10  cm  thick  breasts.  Here,  we  show  3  patients 
who  have  post-operative  lesions  with  optical  images  and  a  corresponding  MRI  (Magnetic  Resonance 
Imaging). 


Case  No.  1 

The  optical  images  shown  in  Figure  2  to  5  are  from  a  38  year-old  patient  Her  cancer  was 
removed  one  month  previously  (lumpectomy),  and  under  pre-  and  post-radiation  therapy, 
respectively.  Figures  2  and  3  show  post-lumpectomy,  pre-radiation  images.  Figmes  4  and  5  show 


Figure  2.  MRI  (lateral  slice)  of  case  1,  whose  cancer  was  removed  (lumpectomy)  one  month 
previously.  This  MRI  was  taken  prior  to  radiation  therapy  (pre-radiation).  Tlie  closed  area  shows  the 
area  where  optical  imaging  was  conducted  (see  Figure  3).  Also  note  that  a  2  x  1.5  cm  seroma  is  seen 
in  the  location  where  the  liraipectomy  had  b^n  perform^ 
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The  lightguide  separation  was  6  cm.  The  seroma  had  less  oxygen  concentration;  O2  saturation  in  Hb 
is  82%  in  the  seroma  and  85-87%  in  the  surrounding  tissue  pre-radiation  (Figure  3-C:).  After 
radiation  therapy,  the  size  of  the  cyst  became  reduced  and  form^  a  heterogeneous  fihrotic  tissue 
(Figure  4),  An  overall  lower  Oi  level  was  observed  in  the  same  area  post-radiation  as  pre-radiation 
(Figure  5C  as  compared  to  Figure  3C),  and  similar  heterogeneity  was  observed  in  the  area  where  the 
smaller  seKMna  is  located.  The  area  where  the  seroma  and  fibrous  scar  tissue  are  located  has  less 
blood  volume  Gess  absorption  (Figures  3-B,5-B))  and  more  scattering  (Figures  3-A,  5- A). 


A/////////////  u'////////////,  c-y////////////. 
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Figure  3.  Optical  images  of  case  1,  under  pre-radiation  corresponding  to  Figure  2.  A:  scattering 
image.  B:  blood  volume  image  C:  oxygen  image.  Note  that  the  seroma  has  low  oxygen 
concentration  than  the  surrounding  tissue. 


Figure  4.  MRI  of  case  1,  after  a  radiation  therapy.  The  closed  area  shows  the  area  that  optical 
imaging  was  conducted  (see  Figure  5^  The  size  of  the  seroma  became  reduced  and  formed  a 
heterogeneous  fibrotic  mass  after  the  irradiation. 
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Figure  5.  Optical  images  of  case  1,  post-radiation  therapy,  corresponding  to  Figure  4.  A;  scattering 
image.  B:  blood  volume  image  C:  oxygen  image.  The  lesion  shows  more  scattering  and  contains 
less  oxygen  than  the  surroun^ng  tissue. 

liie  MRI  and  the  optical  image  (Figures  6,7)  are  from  a  post  lumpectomised  cancer  patient. 
She  is  48  years  old  who  has  small  to  medium  size  breasts.  A  small  lightguide  sep^tion  of  4  cm  was 
used.  After  lumpectomy,  there  is  a  4.5  x  2.5  cm  mass  consisting  of  fibrotic  scar  tissue  and  seroma  in 
the  middle  of  her  breast  (Figure  6).  The  optical  images  were  taken  laterally  covering  between  seipnw 
and  fibrotic  tissues.  The  oxygen  image  shows  that  a  part  of  the  serona  has  a  low  oxygen  level  indi¬ 
cating  a  quite  heterogeneous  profile  figure  7-C).  The  oxygen  level  in  the  seroma  capsule  is  particu¬ 
larly  higher  than  the  surrounding  tissue.  The  blood  volume  image  shows  that  the  seroma  has  less 
blood  volume  (Figure  7-B).  In  addition,  there  is  heterogeneity  in  the  scattering  images  in  the  seroma 
as  well  as  in  the  surrounding  fibrotic  tissue  (Figure  7-A)  constracting  an  image  of  scar  tissue. 


Figure  6.  MRI  of  case  2,  who  has  a  large  lesion  after  a  lumpectomy.  The  lesion  consists  of  seroma 
and  fibrous  scar  tissue. 
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Figure  7.  Optical  images  of  case  2.  The  images  were  taken  in  the  cranio-caudal  direction.  The 
illumination  was  along  the  direction  of  indicated  lines  and  between  two  lines  (Figure  6).  The  line 
nearer  to  the  nipple  corresponds  to  the  left  line  of  each  image.  The  right  side  of  the  images  clearly 
shows  the  lesion.  Note  that  the  lesion  has  a  location  with  lower  oxygen  which  also  scatters  more 
than  surrounding  tissue,  indicating  heterogeneity. 

Case  no.  3 

The  patient  is  52  years  old  and  weighs  180  lbs.  Her  tumor  in  the  breast  has  been  determined 
to  be  cancer  by  biopsy.  She  already  had  been  treated  with  a  lumpectomy  followed  by  radiation 
therapy  (Figure  8).  The  optical  images  (Figure  9)  were  taken  with  cranio-caudal  direction.  In  the 
optical  scattering  image  (Figure  9- A),  there  are  two  highly  scattering  masses,  connected  to  each  other 
by  a  capsule  like  structure.  This  capsule  like  stracture  has  higher  oxygen  tension  than  in  the  two 
masses  (Figure  9-C).  The  two  masses  differ  in  blood  volume  (See  Figure  9-B),  and  the  one  nearer 
the  nipple  fcis  more  blood  volume  than  the  other. 


Figure  8.  MRI  (cranio-caudal  direction)  of  case  3,  post-lumpectomy.  The  closed  area  shows  the  area 
in  which  optical  imaging  was  conducted  (see  Figure  9). 
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Figure  9.  Optical  images  of  case  3,  corresponding  to  Figure  8.  A:  scattering  image.  B:  blood  volume 
image  C:  oxygen  image.  There  are  two  highly  scattering  masses  connected  to  each  other  by  a  capsule 
like  stmcture.  TTiis  capsule  has  higher  oxygen  tension  tfian  in  the  two  masses. 

1.4  Discussion 

Clearly,  we  have  denoonstrated  that  the  time  gated  optical  imaging  technique  is  capable  of 
imaging  a  relatively  large  object.  However,  this  technique  does  not  yet  clearly  detect  a  small  cancCT 
(generally  smaller  than  1  cm),  therefore  we  did  not  include  those  small  tumors  in  this  report.  It  is 
shown  the  timp.  gated  imaging  system  is  not  yet  as  precise  as  X-ray  mammo^phy  for  lesions 
smaiw  than  1  cm  in  size.  However,  it  is  certainly  better  than  the  current  diaphanography  in 
principle. 

In  the  optical  imaging  technique,  the  main  cause  responsible  for  low  resolution  of  an  image  is 
photon  diffusion.  For  the  human  breast  image,  conditions  for  the  continuous  wave  technique  (any 
technology  which  uses  CW  light  ot  current  diaphanography)  aK  viable  in  breasts  that  are  small,  ^ 
and  fatty,  and  having  tumors  larger  than  1  cm  in  diameter.  Since  the  CW  technique  will  not  satisfy 
most  of  the  needs  for  early  detection  of  breast  cancers,  we  have  develop  a  time  gated  photon 
imaging  system  to  test  feasibility  for  the  purpose.  Scattering  and  absorption  images  use  only  15-20% 
and  40-50%  of  early  arriving  photons,  respectively,  which  therefore  eliminates  most  of  the  diffusive 
photons  that  create  obscurity.  The  results  show  that  the  scattering  image  may  have  more  clear  images 
of  lesions  than  those  used  for  absorption  (oxygen  and  blood  volume  images)  because  of  narrower 
time  windows  for  scattering  images. 

It  is  clear  that  this  time  gated  system  is  not  adequate  to  show  a  high  contrast  of  lesions  of  less 
than  1  cm  in  size  from  surrounding  tissue.  The  cause  of  low  resolution  comes  from  mmy  factors. 
First,  although  fortunately  the  breast  contains  substantial  amounts  of  fat,  whose  scattering  factor  is 
low  compared  to  any  other  tissue  constituents,  contrast  between  fat  and  fiber  versus  tunwr  is  not 
always  substantial  enough.  Secondly,  differences  of  oxygen  tensions  or  blood  volumes  between 
those  tissue  components  also  may  not  be  substantial  enough  for  the  available  signal  levels.  Thirdly, 


11 


the  technique  we  use  here  presents  a  number  in  a  pixel  as  a  mean  value  of  a  long  tissue  mass  across 
the  two  lightguides,  and  can  not  represent  a  small  lesion. 

These  problems  can  be  solved  by  finding  causes  of  obscurities.  A  more  sophisticated  image 
reconstructing  technique  such  as  inverse  recovery  technique  (5),  will  better  atl^t  photon  ^fusion 
characteristics  into  the  image  construction  procedures.  In  addition,  we  will  be  able  to  obtain  higher 
resolution  images  with  the  phased  array  system  and  with  a  contrast  agent,  such  as  cardio^een, 
which  may  accumulate  in  the  lesions  widi  permeable  capillaries  and  more  circulation.  We  arc  in  the 
early  stage  of  developing  the  optical  imaging,  and  in  the  future,  attempts  to  improve  substantially  the 
quality  of  optical  imagines  will  be  pursu^ 

2.  TRS  STUDY  OF  OPTICAL  CHARACTERIZATION  (Years  1  and  2) 

2.1  Introduction 

We  have  used  the  TRS  device  to  image  breast  tumors  with  optical  outcomes;  absoiption 
coefficient  and  scattering  coefficient  are  shown  in  Section  1.  At  the  same  time,  the  data  can  be 
analyzed  quantitatively.  The  patient  population  studied  consisted  of  two  major  populations;  tumor 
bearing  breasts  and  post  surgical  breasts  (lumpectomized).  The  purpose  of  separating  lumpw- 
tomized  breast  from  tumor  bearing  breast  was  to  seek  out  some  differences  in  the  post  operative 
breasts  and  tumor  bearing  ones.  It  is  interesting  to  see  if  the  X-Ray-dense  fibrocystic  breast  can  be 
recognized  by  the  optical  means.  In  this  sense  we  separated  the  fibrocystic  disease  with  tumor  frtHn 
simply  tumor  bearing  breasts  even  though  there  were  only  three  patients  in  this  category.  In  the  total 
of  63  patients  studied,  we  are  able  to  analyze  51  patients’  data  at  this  point 

2.2  Methods  (TRS) 

For  the  optical  measurements  using  TRS,  tiie  methods  are  described  in  detail  in  Section  1 . 
The  data  were  taken  usually  at  a  selective  area  of  the  lesion,  2  to  4  cm  square  area  of  the  planes 
created  by  mildly  compressing  side  by  side  or  up  to  down  direction.  TTie  patients  usu^y  had 
between  25  to  64  points  of  measurements  dependent  upon  the  area  of  lesions.  The  opticd  dam 
analysis  are  available  with  fiilly  developed  software.  In  brief,  the  time  resolved  spectrum  is  fitted  in 
diffusion  equation  (6)  with  the  minimum  square  root  method.  The  dam  are  averaged  to  present  only 
with  one  v^ue  from  one  subject 

2.3  Results  (TRS) 

The  breasts  of  51  subjects  were  analyzed  quantitatively  at  this  point  The  patient  populations 
were  divided  into  three  categories;  lumpectomized,  tumor  bearing,  and  fibrocystic  breasts.  Table  I 
shows  the  mean  and  standard  deviation  of  the  pa,  ps’  and  oxygen  saturation  with  hemoglobin  in  the 
normal,  lumpectomized,  tumor  bearing,  and  fibrocystic  breasts.  The  absorption  coefficients  in  die 
normal  population  (n=18)  have  a  lowest  mean  values  (0.025),  and  all  three  patient  populations  have 
double  the  value  of  pa  (see  Figs.  10  and  12).  The  scattering  coefficients  follow  similar  high  v^ues 
to  pa.  In  the  fibrocystic  breast,  ps’  was  much  greater  than  any  other  group,  i.e.,  three  time  higher 
than  the  normal  breasts  in  the  three  subjects  we  measured  (see  Figs.  1 1  and  13). 
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Fig.  10.  Histogram  of  the  absorption  coefficient 
with  all  patient  populations  and  normal  controls. 
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Fig.  11.  Histogram  of  scattering  with  all  patient 
populations  and  normal  controls. 
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Fig.  12.  Histogram  of  p,:  absorption  coefficient 
with  patients:  tumor  bearing  and  normal  control. 
A  greater  absorption  coefficient  is  seen  in  the 
patient  population. 
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Fig.  13.  Histogram  of  p,’  in  the  patient  popula¬ 
tions:  tumor  bearing  breast,  and  normal 
controls. 


The  above  pa  histogram  is  a  display  of  the  tumor  bearing  breast  (cancCT  and  benign  tumor)  as 
compared  to  normal  controls.  A  greater  pa  is  seen  in  the  tumor  bearing  populations.  The  histogram 
of  ps’  is  of  the  tumor-bearing  breasts  (cancer  and  benign)  as  compared  to  normal  control  subjects. 


2.4  Discussion 

We  expected  that  because  tumor  growth  is  usually  accompaiued  by  angiogenesis  and  higher 
blood  flow,  high  blood  content  could  be  seen  and  it  seems  that  is  the  case  presented  in  this  study. 
Also,  lumpectomized  breasts  within  2  weeks  should  still  produce  higher  metabolism  because  of  tiie 
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healing  process,  leading  to  such  high  absorption  coefficients  due  to  high  blood  contents.  It  will  be 
interesting  to  see  the  effect  of  time  after  surgery  on  pa. 

It  would  be  \ery  interesting  if  we  can  analyze  more  precisely  those  pa  tmd  ps’  values  related 
to  the  specific  lesion.  However,  we  do  not  know  the  exact  location  of  the  lesion  at  the  time  of  the 
optical  measurement  which  prevents  the  tumor  from  being  analyzed  with  regard  to  pa  and  ps’  more 
closely.  Therefore,  only  one  way  to  analyze  the  tumor  characterization  is  during  another 
measurement  such  as  MRI  and  X-Ray  mammography.  That  is  exactly  what  we  did  in  the  trial  for  the 
last  phase  of  this  study.  It  is  not  surprising  to  see  Ae  hyperemia  near  the  tumor  tissue,  within  the  2.5 
array.  Also,  it  may  be  trae  that  in  a  small  tumor,  this  hyperemic  area  may  be  much  larger  than  the 
area  of  tumor  itself. 

Further  analysis  of  the  rest  of  unanalyzed  patients  and  controls  may  provide  statistically 
significant  results  for  pa  and  ps’.  The  increase  in  pa  can  indicate  some  unusual  breast  metabolic 
activity,  which  can  be  used  for  screening  purposes. 

3.  PRESENT  AND  FUTURE  PLANS  FOR  THE  TRS  SYSTEM  FOR  OPTICAL 

CHARACTERIZATION  IN  TUMOR  BEARING  BREASTS 

The  detailed  tasks  undertaken  in  1995-6  are  in  accord  with  the  1994  proposal:  the 
development  of  technology  for  in-magnet  recording  of  photon  migration  data  on  human  breast  is 
considered  mandatory  for  the  following  reasons: 

1.  Co-registration  of  MRI  and  TRS  imaging  is  immediately  possible. 

2.  Access  to  the  patient  population  and  to  facilitated  patient  consent  is  assured. 

3.  A  priori  knowledge  of  the  location  of  the  tumor  from  (and  X-ray  mpimographic) 
coordinates  facilitates  the  imaging  of  the  tumor  and  the  calculation  of  its  scattering  (ps')  and 

absorption  (pa)  properties. 

4.  Advanced  algorithms  have  now  been  developed  for  creation  of  images  of  tumora  from 
several  milHmeters  to  bettCT  than  1  cm  in  resolution  and  as  mentioned  above  for  cal^ating  tiie 
tumor’s  optical  properties.  Such  algorithms  are  described  in  the  Ph.D.  Theses  of  Dr.  David  Boas  (6) 
and  Dr.  Maureen  O’Leary  (7)  and  in  a  recent  contribution  by  Xingde  Li  (8).  All  these  algorithms  can 
effectively  be  applied  to  the  experimental  data  gathered  in  the  NMR  magnet. 

5.  Sensitivities  available  to  the  optical  mefliod:  The  availaWlity  of  pa  and  ps'  data  for  the 
tumor  volume  and  for  the  adjacent  area  enables  the  calculation  of: 

a.  Increased  blood  concentration  in  a  tumor  due  to  hypervascularization. 

b.  Decreased  oxygen  saturation  in  the  tumor  due  to  hypermetabolic  state  of  the  tumor  with 
respect  to  the  surrounding  breast  tissue. 

c.  Intensity  response  to  cardiogreen  (due  to  indocyanine  green,  ICG).  This  is  a  vascular 
probe  which  not  only  indicates  imnoediately  the  vascular  volume,  but  ovct  time,  indicates  the  extent 
and  course  of  its  extravasation,  generally  used  now  with  MRI  to  indicate  malignancy. 

Light  scattering  and  absorption  of  tumor  bearing  human  breast  has  been  studied  in  a  large 
number  of  X-ray  mammographicaUy  studied  patients,  'nme  domain  data  acquired  using  the 
methodology  of  Fig.  14  enabled  a  fi^quency  domain  analysis  to  be  made  by  Foinier  transformation. 
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3.1  Experimental  Methods/Results  Obtained 

Transformation  from  Time  to  Freauencv  Domain.  In  the  foregoing  tasks  we  have  propos^ 
both  time  domain  and  frequency  domain  methods  to  be  studied  and  to  be  evaluated  for  the  effic^y  in 
breast  tumor  detection  (9).  Furthermore,  these  methods  have  required  significant  development  in  die 
course  of  the  research  and  finalization  of  the  designs  of  both  these  methods  has  been  achieved. 
Preliminary  data  has  been  obtained  on  model  systems  which  can  afford  the  basis  for  extending  the 
studies  to  an  appropriate  size  of  tumor  bearing  population  as  had  been  contemplated  in  Table  I  of  the 
original  application  and  where  in  the  year  1  repoii  the  three  methods  were  described  ((See  Appendix 
I  for  Table)  and  a  series  of  patients  expected  to  be  studied.  In  order  to  satisfy  the  needs  of  this  table, 
over  1,500  studies  were  made.  A  special  report  was  made  of  23  observations  on  the  oxygen 
saturation  in  the  tumor  bearing  breast  (see  below).  Furthermore,  histogram  displays  of  pa  and  ps’  for 
human  breasts  with  small  tumors  wCTe  made  to  establish  the  distribution  of  values  to  be  expected 
using  the  “banana”  pattern  optical  detection  of  tumors. 
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Fig.  14.  A  detailed  diagram  of  the 
components  for  in-magnet  TRS  Optical 
Tomography  system.  The  clock  pulses 
trigger  the  laser  and  the  TAC.  The  time 
selector  affords  variable  timing  of  the 
selection  of  fibers  for  general  search 
and  for  localized  quantitation.  Parallel 
detection  of  the  optical  signals  together 
with  parallel  processing  achieves 
optimal  signal  to  noise  ratio.  Photon 
ouqtuts  through  the  breast  are  adjusted 
in  intensity  for  each  position  of  the  light 
source  to  ensure  optimal  photon 
counting  conditions.  Computers  for 
image  reconstraction  are  provided. 
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Year  2  reported  detailed  studies  of  a  total  of  three  groups  of  patients,  which  gives  a  tot^  of 
140  examinations  v^dating  the  required  studies  of  Table  1  (Appendix  I).  23  subjects  were  examined 
for  oxygen  saturation  which  showed  that  50%  of  the  breasts  to  be  in  a  hypoxic  state,  and  therefore, 
potenti^ly  radiation  resistant 

Future  plans  for  Year  2  involved  the  setting  up  simultaneous  galdolium  enhanced  and 
cardiogreen  enhanced  optical  tomography  in  order  to  obtain  greater  specificity/sensitivity,  and 
indeed,  co-registration  of  the  images.  ITiis  has  necessitated  substantial  instrumentation  developinent 
according  to  Fig.  2  of  the  Year  2  Report  (Appendix  11)  for  which  additional  funds  wwe  recruited 
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(Fig.  14  above).  Furthermore,  a  phase  modulation  system  was  developed  in  ordw  that  1  to  1 
comparison  of  the  two  systems  might  be  available  for  further  work  under  Army  Contract  support. 

3.2  Co-Registration  System:  MRI/Optical  Imaging 

The  preliminary  diagram  of  this  system  presented  in  Fig.  2  of  the  Year  2  (Appendix  D)  report 
is  now  evolved  as  in  Fig.  14  above.  The  system  consists  of  a  16-channel  source  detector  system 
with  a  multiwavelength  laser  diode  illumination.  This  illumination  is  time  sequenced  through  16 
fibers  to  successively  illuminate  the  breast  in  16  positions  at  multiple  wavelengths,  usually  760  and 
780  nm  and,  in  some  cases,  830  nm.  The  laser  diode  power  is  wito  the  FDA  limits  of 
approximately  20  microwatts  and  is  delivered  in  pulses  of  20  to  100  pico  sectMid  duration.  16 
detectors  simultaneously  observe  the  photon  migration  signals  by  transmission  through  the  b^t 
and  measure  the  photon  kinetics  at  the  various  positions  at  which  available  signal  to  noise  ratio  is 
obtainable.  The  16  channels  are  selectively  processed  by  a  multiplexer  router  (Edinburgh  Electronics 
or  EEG  Electronics)  and  thus  multichannel  analysis  is  available  at  the  time  course  at  the  selected 
wavelengAs  and  at  the  selected  positions.  The  data  on  photon  migration  is  stored  in  a  computer 
which  may  be  analyzed  according  to  the  algorithms  used  previously  in  the  year  2  survey,  or  by 
imaging  algorithms  which  have  been  developed  by  Boas  and  O’Leary  (6,7)  and  are  discussed  as  a 
part  of  this  final  report 

Theimplementationof  this  device  in  the  1.5  T  magnet  is  illustrated  in  Fig.  15  which  shows 
the  orientation  of  the  16  source/detector  combinations  wiA  respect  to  the  NMR  coils.  Furthermore, 
small  capillaries  are  filled  with  water  and  1  millimolar  copper  sulfate  in  order  to  mark  clearly  die 
relationship  between  the  NMR  coils  and  the  optical  fibers. 


Fig.  15.  Optical  coupling 
to  the  soft  compression 
plates  for  MRI  imaging 
as  shown  by  the  open 
circles  for  sources  and 
the  solid  circles  for 
detectors.  The  fibers  will 
lie  parallel  to  the  plates 
because  of  the  right 
angled  prisms. 


3.3  Current  State 

This  apparatus  has  been  assembled  using  8  channels  and  a  microchannel  plate  detector.  16 
PMT’s  (photomultiplier  tubes  T08)  have  been  ordered  fi-om  Hamamatsu  and  will  be  delivered 
shortly.  This  is  more  efficient  than  the  microchannel  plate  and  has  less  crosstalk. 
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4.  CWS  BREAST  IMAGING  WITH  CONTRAST  AGENT  (Years  2  and  3) 

4.1  CWS  Methods 

Eight  or  16  combinations  of  li^t  sources  (tungsten  bulbs)  and  light  detectors  (silicon 
photodiodes)  are  located  on  a  two  dimensional  circle  line  equal  distance  apart  from  each  other  (22.5° 
or  45°  angle  apart,  Figures  16A,B). 

~  from  source 

1 


to  detectors 

Figure  16.  The  scheme  A  describes  the  cut  section  of  the  16  light  source-detector  device,  holding  a 
human  breast  inside.  The  diameter  can  be  fitted  easily  with  set  screws.  When  the  device  fits  a  breast, 
source  detector  relation  remains  unchanged  (B). 

4.2  Image  Algorithm,  Back  Projection 

The  original  back  projection  algorithm  has  been  developed  for  the  CAT  scanned  data.  Since  that 
application  deals  with  straight  rays,  when  it  is  applied  to  diffusive  rays  of  photons,  inodifications 
must  be  to  adapt  the  characteristics  of  the  photons.  Firstly,  the  idea  of  "photon  hitting  density" 
(5)  replaced  the  linear  relationship  of  X-ray  when  the  locations  of  be^  are  trpslated  into  pixels. 
Photon  hitting  density  is  the  probability  of  a  photon  occupying  a  specific  place  in  between  a  source 
and  a  detector.  Figure  17  shows  the  area  where  photons  can  be  localized  with  approximately  80  to  60 
%  probabilities  when  it  is  assumed  that  pa  and  its'  are  known  to  be  similar  to  the  human  breast  and 
homogeneous. 


Figure  17.  Photon  hitting  density  image  of  breast.  Approximately  80(A)  to  60(B)  %  of  the  photons 
will  pass  through  these  shaded  area. 

In  the  image  re-construction  program,  the  probability  is  translated  into  a  weight  factor  when  it  is  used 
to  process  back  projection.  The  back-projection  is  basically  a  process  averaging  out  values  of 
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information  that  each  beam  carries  with  the  weighing  in  each  pixel.  Therefore  we  can  not  recovw 
correct  values,  which  is  attenuated  or  amplified  by  the  surrountfings. 

4.3  CWS  Model  Study 

In  the  circle,  we  placed  a  breast  model  using  0.5%  intralipid  and  ink  (about  the  same  scattering  and 
absorption  coefficients  as  those  in  the  breast  tissue).  It  takes  only  a  few  seconds  to  acquire  data  since 
one  source  is  on.  All  the  detectors  detect  the  photons  at  the  same  time.  16  soirees  cycle  at  a  rate  of 
200  msec/source  (Fig.  18). 


Figure  18.  The  model  study  of  optical  imaging  with  the  8  probes.  Intralipid  0.5%  and  ink  resembles 
the  optical  characteristics  of  the  human  breasts.  Two  1  cm  objects  are  located  1  cm  apart.  When  8x8 
images  are  made,  the  two  objects  are  depicted  in  the  reconstructed  image  by  the  back  projection 
algorithm. 


The  data  acquisition  for  an  image  takes  4-8  seconds  with  8x8  or  16x16  data  points.  A  contrast 
agent;  Indocyanine  Green  (ICXj,  0.25  -  0.5  mg/kg),  is  injected  in  the  model  system  similar  to  treating 
human,  when  data  are  being  acquired  continuously  every  few  seconds  over  10  -  15  min  as  the  ICG 
absorption  signal  appears  and  disappears.  Data  are  processed  as  OD  changes  due  to  ICG  from  the  iv 
injection.  Some  data  points  are  averaged  in  several  time  windows  and  as  a  result,  several  images  are 
reconstructed  using  back-projection  algorithm.  The  results  show  that  with  the  ICG  concentration 
which  we  inject  into  humans,  apparent  resolution  is  7  mm  for  the  contrast  agent 

4.4  Human  CWS  Imaging  with  Contrast  Agent 

The  patient  (56  years  old)  has  a  long  history  of  a  palpable  tumor.  The  lump  swells  with  menstmation 
cycle.  When  ICG  is  injected,  optical  image  depicts  the  lesion  area  (Fig.  19,  MRI),  and  contrast  of  the 
lesion  are  great  at  the  following  minutes  (2-4  minutes.  Figs.  20C-6E)  in  the  16x16  image.  It  reveals 
that  a  functional  tumor  like  tiiis  also  has  a  greater  circulation  and  is  associated  with  a  high 
accumulation  of  ICG.  The  contrast  agent  makes  it  possible  to  depict  the  tumor  in  this  case. 

In  the  same  patient,  MRI  shows  a  colony  of  small  fibrous  cysts,  with  ducts  which  are 
hypertrophic  and  connected  to  the  nipple  (see  Fig.  20  MRI). 
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The  results  show  diat  the  resolution  of  tumor  detection  is  enhanced  a  great  deal  (note  the 
comparison  between  Figs.  19A  and  19C).  It  also  shows  that  wme  of  the  benign  tumor  will  be 
detected  and  further  diagnosis  has  to  be  made  after  this  type  of  image.  It  is  possible  that  we  could 
find  a  way  to  distinguish  a  malignant  tumor  by  the  optical  characteristics,  such  as  by  a  combination 
of  pa  and  ps'  in  malignant  tumor,  or  contrast  agent  time  course  kinetics  due  to  leakage  of  vessels 
with  high  angiogenesis  activity  etc.  These  are  still  yet  to  be  studied. 
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Figure  19.  Optical  images  of  a  breast  with  fibrocystic  disease  with  ICG  injection  at  B. 

A  -  pre-injection.  B-F  correspond  to  one  minute  images  after  the  ICG  injection  (1-5  min). 


Figure  20.  MRI  of  a  breast  with  fibrocystic  disease.  Gadolinium  accumulation  is  seen  in  the  swelling 
ducts  of  the  cyst. 


5  IMAGING  USING  THE  FREQUENCY  DOMAIN  METHOD  (Years  2  and  3) 

5.1  Phased  Array  Imaging 

The  rapidly  modulated  light  obtained  from  a  laser  diode  can  be  employed  equally  well  and 
affords  a  somewhat  more  contact  system  for  breast  tumor  detection.  In  this  case,  the  breast  is 
immersed  in  matehing  intralipid  as  shown  in  Fig.  21  and  a  dual  wavelength  scanner  is  eleclro- 
mechanically  passed  by  the  breast  immersed  in  the  above  mentioned  matching  fluid.  The  tank  is  wide 
enough  to  accommodate  the  larger  than  average  breast. 


Fig.  21.  Illustrating  a  breast  scanner  with  pendent  breast  (in  Saran  wrap)  immersed  in  matching 
intralipid,  ink  p,  =  0.02,  p,’  10  cm*‘  (warmed  to  body  temperature,  sources  and  detector  translated 
along  box  at  appropriate  heights). 

The  performance  of  this  device  is  indicated  in  scanning  of  model  breast  which  is  occupied  by 
the  equivalent  of  a  half  milligram  per  kilo  cardiogreen  (ICG)  compound.  The  ICG  is  the  choice  of 
agent  for  use  in  our  studies  as  a  complement  to  the  galdolium  chelate,  and  is  also  useful  for  the  MRI. 
In  this  case,  the  trace  shows  the  detection  of  the  location  of  a  1  cm  diameter  tube  containing  6 
miciomolarICG  so  that  the  active  volume  seen  by  the  prote  is  approximately  12  nanomoles  of  die 
cardiogreen,  an  amount  that  can  be  expected  to  be  obtained  in  a  1  gram  tumor  (see  Fig.  22). 


Fig.  22.  Illustrating  the  accuracy  of  location  of  a  1  cm^  object  filled  with  6  pM  ICG.  The  object  is 
located  at  the  zero  volts  intersection  +  1  mm. 
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A  PMS  system  (50  MHz,  NIM  Inc.)  with  5  mwatt  780  nm  laser  diodes,  having  two  phase  O’ 
and  180°  angle  was  used  for  imaging.  These  two  lasers  are  located  on  the  circular  surface  at  22°  or 
45°  angle,  while  the  detecting  lights  from  360°  angles  with  lightguides  16  or  8  respectively,  similar  to 
the  CWS  imager  (Figure  23).  A  tumor  model  with  intralipid  0.5%  and  black  plastic  rod  cylinder  was 
placed  in  the  intrepid  container. 


Figure  23.  PMS  image  of  single  object  model  (black  plastic  rod  0.8  cm)  with  0.5%  intralipid 
surrounding  in  a  8  cm  container.  (Left)  Schema  for  8  probe  imager.  Two  sources  (SI,  S2)  are 
located  at  45°  from  each  other.  Detector  6  receives  phase  shift  from  the  two  light  sources.  A.  image 
using  8x8  imager.  B.  image  using  16  x  16  imager. 

5.2  Imaging  Developments:  FFT  Approach  to  Biomedical  Imaging  with  Diffusive 
Photon  Density  Waves. 

In  recent  years,  a  variety  of  techniques  for  imaging  turbid  media  such  as  tissue  with  diffuse 
light  have  been  explored.  Most  of  these  methods  employ  direct  matrix  iriversion  (e.g.,  singular  value 
decomposition)  or  iterative  techniques  (e.g.,  ART.  SIRT  and  Conjugate  Gradient)  for  image 
reconstracticMi,  and  they  are  in  genei^  computaticMially  intensive.  In  this  paper,  we  introduce  a  new 
imaging  methc^ology  which  is  computationally  much  faster.  Our  ^proach  is  essjwitially  a  near-field 
wave  technique  that  relies  on  a  series  of  2-dimensional  fast  Fourier  transforms  (FFT’s).  We  present 
a  rigorous  account  of  Ae  theory,  and  the  first  experimental  images  of  absorbing  and  seating 
objects  in  turbid  media  using  the  approach.  In  addition  to  producing  information  about  the  position 
and  shape  of  the  hidden  object(s),  we  have  found  that  under  some  circimstances  it  is  possible  to 
ntiliTf.  projection  images  to  deduce  the  optical  properties  within  a  thin  slice  without  the  need  for  a 
complex  reconstmetion  procedure  such  as  matrix  inversion.  It  should  be  possible  to  obtain  clinical 
projection  images  in  real  time  widi  this  fast  FFT  approach.  (10)” 

This  method  (See  Appendix  HI)  has  recently  been  used  to  obtain  images  of  two  normal  adult 
human  breasts,  in  one  case  wiA  the  tumor  phantom  in  the  breast  chamber  of  Fig.  21,  as  presented  at 
the  SPIE  1997  meeting. 

In  this  study  the  cancer  was  not  seen  in  the  X-ray  mammography,  but  was  palpable  (^ 
size).  The  image  was  reconstructed  primarily  by  a  near-field  wave  technique  that  relies  upon  a  series 
of  2-dimensional  FFT’s  as  described  above. 
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Figure  24.  The  breast  image  presents  pa  (calculated  from  the  frequency  domain  system)  on  a  patient 
with  a  cancer  close  to  her  chest  (upper  right) 


CONCLUSION 

We  demonstrate  here  that  with  three  different  technologies,  we  c^  de^  cancers.  It  is 
however  premature  to  conclude  resolution  limitations  of  each  technique  at  this  point.  The  resolution 
we  presented  here  in  TRS  images  is  in  a  range  of  1  cm,  but  this  photographic  image  is  not  a  valid 
way  to  judge  resolution,  because  the  thickness  of  the  breast  has  not  been  resolved  in  the  current 
method.  However,  this  study  has  shown  that  it  should  be  much  better  than  1  cm  resolution  when  die 
thickness  of  the  breast  can  te  resolved.  Surprisingly,  the  diffusive  CWS  has  a  potential  for  imaging 
human  breast  lesions.  The  key  for  the  success  is  to  produce  higher  contrasts  by  the  use  of  contrast 
agents  with  the  modified  back-projection  algorithm.  However,  we  show  here  an  example  of  die 
benign  tumor,  which  clearly  indicates  greater  blood  flow  and  blood  volume  in  the  hyperactive  benign 
tumor.  A  serious  question  here  is  does  the  blood  flow-  bloodA'olume  image  really  diagnose  a 
malignant  tumor?  We  Can  only  answer  this  question  by  studying  more  cancer  cases  as  well  as  benign 
tumors  and  compare  these  two. 

The  next  step  is  to  search  for  the  specificity  of  a  malignant  tumor.  We  now  know  that  to 
image  the  high  circulation,  blood  volume  does  not  necessarily  identify  the  malignancy  because  it 
shows  a  hyperactive  metaltolism  or  the  heterogeneity  of  the  tissue  metabolism.  The  other  function  of 
high  circulation  high  blood  volume  image  has  to  te  considered.  For  example,  an  unusually  high 
permeability  constant  or  diffusion  coefficient  of  ICG  into  the  lesion  may  identify  n^gn^cy.  'Hie 
contrast  method  we  used  may  help  identify  malignancy  by  observing  the  time-distribution 
relationship  of  the  contrast  agent 


22 


Many  previous  investigation  have  indicated  that  oxygen  concentration  in  cancer  is  tnwh 
lower  than  in  the  normal  tissues  (11).  Benign  tumors  do  not  have  such  low  oxygen  concentration 
Ul).  The  oxygen  image  should  be  considered  as  a  potential  diagnostic  image  for  cancer.  Oxygen 
images  can  be  easily  obtained  by  adding  the  information  of  multiple  wavelength  data. 
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Appendix  I 

Table  from  Original  Application  and  Year  1  Report 


PLANNED  SCHEDULE  OF  STUDIES 


Year 

1 

2 

3 

Methods 
(Pulse  Time) 

+ 

- 

- 

Methods 

(Phase  Modulation) 

- 

+ 

+ 

Method  4* 

(Amplitude  Modulation) 

- 

+ 

+ 

No.  of  Studies 

120 

120 

120 

Expected 

Malignant 

so 

so 

so 

*120  patient  studies  with  this  method  are  proposed.  It  is  proposed  that  50  of  these  be 
studi^  in  Years  2  and  3. 


Appendix  n  -  From  Year  2  report 


Future  Plans:  Aims  2^,  and  4. 

As  is  consistent  with  Aims  2, 3  and  4,  we  are  proposing  as  a  part  of  the  Year  3  effort  md  as  a 
main  feature  of  any  possible  further  funding  to  carry  out  the  optical  and  MRI  examin^oiis 
simultaneously  rather  than  consecutively  since  only  a  small  fiactiCMi  of  the  available  population  is 
willing  remain  to  continue  the  study  with  the  optical  method.  The  diagram  of  Figure  2  illusirates  our 
proposal  to  accompUsh  this  by  fiber  optical  coupling  of  the  optical  method  direcdy  into  the  magnet 
bore  through  34  feet  of  32  fibers,  16  on  each  side  of  the  breast,  and  an  additional  16  fibers  on  the 
contralateral  breast  as  soon  as  the  MRI  has  completed  the  “dual  breast”  coil  assembly,  leading  to  a 
potential  for  64  optical  fibers  for  coupling  in  and  out  of  the  MRI.  Enthusiastic  support  of  this  project 
is  evinced  by  the  MRI  t«»am,  particularly  Dr.  Schnall  who  has  indicated  that  the  soft  compression 
plates  used  with  the  MRI  are  optically  transparent  and  could  readily  be  used  to  support  at  least  16 
source  detector  fibers  on  each  side  of  a  single  breast,  or  in  fact,  both  breasts,  beca,use  MRI  plans  to 
observe  both  breasts  simultaneously  in  order  to  ensure  that  no  tumors  are  located  in  the  connlateial 
breast.  Thus,  sufficient  fibers  will  be  coupled  directly  into  the  MRI  soft  compression  plates  to  afford 
tumor  detection  over  the  entire  volume  imaged  by  MRI,  and  thus  ensure  correlation  of  the  results  of 
the  two  methods.  Furthermore,  the  two  methods  are  non-interfering,  and  optical  data  can  be  acquired 
simultaneously  with  MRI  data  since  no  detectable  magnetic  field  is  induced  by  propagation  of  light 
through  tile  optic  fibers. 


Output  to  other  Hospital  Terminals 
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Figure  2.  Illustrating  a  method  of  procedure  for  fiber  optic  coupling  of  the  time  domain  and  phase 
modualtion  optical  method  to  the  plastic  breast  ctrapression  plate  used  in  the  MRI  studies.  Die 
diagram  indicates  the  use  of  different  fibers  in  the  fiber  optic  coupling  to  obtain  data  from  different 
places  on  the  breast. 
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Abstract 

The  spatial  structure  of  optically  heterogeneous  turbid  media  is  probed  with 
diffusive  hght.  Projection  images  are  obtained  experimentally  by  deconvolut- 
ing  the  scattered  diffuse  photon  density  waves  on  a  planar  boundary  using 
a  Fast  Fourier  Transform  (FFT).  The  method  is  very  fast,  permitting  ob¬ 
ject  locahzation  and  characterization  in  ~1000  volume  element  samples  on 
sub-second  computational  time  scales.  The  optical  properties  of  shce  shape 
inhomogeneities  are  accurately  determined. 

In  recent  years  the  application  of  near-infrared  (NIR)  diffusing  light  for  biomedical  diag¬ 
nosis  and  imaging  has  gained  favor  as  a  result  of  its  non-invasive  nature,  economy,  and  novel 
contrast  relative  to  other  diagnostics  [1-3].  To  this  end  a  variety  of  techniques  for  imaging 
with  diffuse  light  have  been  explored  [4,5].  Most  of  these  methods  employ  direct  matrix 
inversion  (©'g-,  singular  value  decomposition)  or  iterative  techniques  (e.g.,  ART/SIRT  [6] 
and  Conjugate  Gradient  Descent)  for  image  reconstruction. 

In  this  paper  we  introduce  a  new  near-field,  diffusive  wave  imaging  methodology  using 
techniques  similar  to  those  of  conventional  diffraction  tomography  [/].  While  the  near  field 
diffraction  tomography  method  has  attracted  the  attention  of  a  few  researchers  in  the  pho- 
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ton  migration  field  [8],  our  paper  presents  a  rigorous  account  of  the  theory,  and  provides 
the  first  experimental  images  of  absorbing  and  scattering  objects  in  turbid  media  using  the 
approach.  It  differs  from  least-square  techniques  such  as  ART  and  SIRT  in  that  it  is  fast  and 
non-iterative.  In  addition  to  providing  information  about  the  position  and  shape  of  hidden 
object(s),  the  projection  images  can  be  used  to  deduce  the  optical  properties  of  hetero¬ 
geneities  without  the  need  for  complex  reconstruction  procedures  such  as  matrix  inversion, 
when  the  heterogeneities  are  thin  and  information  about  their  depth  is  available. 

We  adopt  the  frequency  domain  picture  for  our  discussion.  An  intensity  sinusoidally 
modulated  light  source  coupled  into  a  highly  scattering  medium  such  as  tissue,  produces  a 
diffuse  photon  density  wave  (DPDW)  that  propagates  outwards  from  the  source  [9].  The 
amplitude  and  phase  of  this  DPDW  depends  on  the  absorption  and  scattering  coefficients 
within  the  turbid  medium. 

In  uniform  turbid  media,  the  DPDW  from  a  point  source  at  Tg  detected  at  position  ru 
has  the  form  t/o(rd,  rs,a;)  =  Mq  exp(i  kolr  -  rs|)/(47rZ)olr  -  rs|).  Here  io  is  the  angular 
source  modulation  frequency,  Mo  is  the  ac  amplitude  of  the  source,  Do  —  u/(3//'o)  is  the 
diffusion  coefficient  {v  is  the  speed  of  light  in  the  medium),  ko  =  yJ{-HaoV  4-  icu)/Z)o  is  the 
DPDW  wave  number,  and  fiao  and  are  respectively  the  homogeneous  absorption  and 
reduced  scattering  coefficients  of  the  medium.  In  heterogeneous  media,  the  total  DPDW, 
t/((rd,  rs,u;),  is  a  superposition  of  incident  (  f7o(rd,rsitu)  )  and  scattered  (  f7i(rd,  rs,u;)  ) 
DPDW’s.  To  the  first  order  in  the  variation  of  optical  absorption  and  reduced  scattering 
coefficients,  the  scattered  wave  is 

[/i(rd,rs,u;)  =  j T{r,Uo{r,rs,^v))  G{\rd  -  r\)  cPr  .  (1) 

V 

Here  (?(|rd  —  r|)  =  exp{fA:o|r  —  rsi)/(47r|r  —  Tgl)  is  the  Green’s  function  for  the 

DPDW  in  the  homogeneous  medium.  We  call  T{r,Uo{r,rs,io))  the  inhomogeneity 
function.  ra6s(r,  f/o(r,  rg,t<;))  =  {—6[ia{r)  v/Dq)  C/o(r,  rs,a;)  for  absorbing  objects,  and 
T,catt{rMo{r,rs,oj))  =  2>Do  ko^ Iv)Uo{v,Ys,uj)  -  V  ln{6n'^  +  •  Vi7o(r,rg,u;)  for 

scattering  objects.  The  integral  is  over  the  sample  volume  V. 


2 


In  this  paper,  we  consider  a  parallel-plane  geonnetry  as  shown  in  Fig.  1(a)  (which  is 
potentially  applicable  to  the  compressed  breast  geometry).  For  this  case,  a  natural  basis 
set  for  the  Green’s  function  in  Eq.(l)  is  the  simple  Weyl  expansion  form  in  terms  of  spatial 
frequencies  p,  m  [10],  i.e. 

z  f  f  dx)  do 

G(|rd  -  r|)  =  ^  J  J  —^exp[ip{xd  -  x)  -|-  iq{yd  -  y)  +  im{zd  -  z)]  ,  (2) 

— oo 

where  we  assume  Zd  >  z  without  losing  generality,  m  =  {k^  —  —  q'^)^  with  Im(m)  >  0. 

Using  Eq.(2)  and  taking  the  2-D  spatial  Fourier  transform  of  both  sides  of  Eq.(l)  (with 
respect  to  transverse  (x,  y)  coordinates),  we  obtain  the  following  relation 

2d 

Ui{p,q,Zd,r^,<^)  =  ^j  f  {p,q,z,  rs,uj)  exp[im{zd  -  z)]  dz  .  (3) 

The  left  hand  side  of  Eq.(3)  is  the  2-D  Fourier  transform  of  the  scattered  DPDW  measured 
on  the  detection  plane  2  =  Zd-  This  integral  equation  is  approximated  by  a  summation 

N 

T{p,q,Zj,rs,io)  exp{-imzj)  = 

where  6z  is  the  discretized  step  size  and  N  is  the  total  number  of  slices  in  the  z  direction. 

For  the  projection  image,  we  replace  Zj  on  the  left  hand  side  of  Eq.(4)  with  the  estimated 
slice  position  of  the  object.  We  drop  the  sum  over  all  other  Zj's  and  then  perform  a  2-D 
inverse  Fourier  transform  of  T  to  obtain  the  projection  image.  When  the  object  thickness  is 
of  order  of  several  transport  mean  free  paths  {1  / {pa  +  p's))  we  are  able  to  deduce  accurately 
object(s)  optical  properties.  For  thicker  objects  (i.e.  >5  mm),  the  average  over  the  size 
of  the  object  weighted  by  the  sum  of  exponential  phase  factors  reduces  the  accuracy  of 
the  optical  properties.  However  position  information  is  accurate,  and  the  relative  optical 
properties  of  multiple  objects  are  also  very  accurate. 

To  demonstrate  the  feasibility  of  this  algorithm,  we  have  performed  amplitude  and  phase 
measurements  in  a  parallel-plane  geometry  (Fig. 1(a)).  The  experimental  setup  is  shown  in 
Fig. 1(b).  The  system  consists  of  an  RF  modulated  (100  MHz),  low  power  (~3  mW)  diode 
laser  operating  at  786  nm.  The  source  light  is  fiber  guided  into  a  large  fish  tank  of  50L 
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i  Sz 


exp{-imZd)Ui{p.  q,  Zd,  Ts, co) 


(4) 
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0.75%  Intralipid  (^^=0.03  cm~\  cm~^)  enabling  the  nse  of  infinite  medium  boundary 

conditions.  A  detection  fiber  couples  the  detected  diffusive  wave  to  a  fast  avalanche  photo¬ 
diode  (APD).  A  single-side-band  (SSB)  demodulator  is  used  to  homodyne  the  signal  and 
the  reference  wave  at  100  MHz.  The  dynamic  range  of  our  current  setup  is  ~2500.  The 
source  and  detection  fiber  optics  are  moved  by  automated  stepper  motors.  The  system  is 
very  stable  for  3  hours. 

The  experimental  geometry  is  shown  in  Fig.  1(a).  The  source  position  was  fixed  and 
taken  to  be  the  origin  of  our  coordinate  system.  As  shown  in  Fig.  1(a),  we  “made”  the 
detection  plane  by  scanning  a  single  detection  fiber  over  a  square  region  from  (-4.65,  -4.65, 
5.0)  cm  to  (4.65,  4.65,  5.0)  cm  in  a  plane  at  Zd=5.0  cm  in  steps  of  size  Ax=Ay=0.3  cm. 
The  amplitude  and  phase  of  the  DPDW  were  recorded  at  each  position  for  a  total  of  1024 
points.  We  directly  measured  the  amplitude  and  phase  in  the  homogeneous  medium  to 
obtain  t/o(rd,  i*s,a;). 

Two  absorbing  slices  each  of  dimension  1.5  x  1.5  x  0.4  cm^  were  then  submerged  in  the 
turbid  medium  (0.75%  Intralipid).  The  slices  were  made  of  resin  plus  Ti02  and  absorbing 
dye.  Slice  1  with  ^ai  =  0.20  cm“^  was  placed  at  position  (-1.6,  -0.3,  3.0)  cm  and  slice  2 
with  iJ.a2  =  0.10  cm""^  was  placed  at  (1.6,  0.3,  3.0)  cm.  The  scattering  coefficients  of  these 
two  slices  are  the  same  as  background,  e.g.,  8.0  cm“h  The  scattered  wave  f/i(rd,  rs,u.') 
was  obtained  by  subtracting  the  homogeneous  background  DPDW  f  oil’d, fs,*^)  from  the 
measured  signal  I7t(rd,  rs,u;).  The  2-D  Fourier  transform  of  7/i(rd,  Ts,^')  leads  to  the  inho¬ 
mogeneity  function  T{p,q,  Zj,rs,uj)  in  Eq.(4).  In  this  paper,  we  use  a  priori  information 
about  the  object  position(s)  in  z-direction  to  select  a  single  image  slice,  e.g.,  a  slice  at  z—Zobj 
where  the  inhomogeneity  function  is  T{p,  q,  Zobj,  fs,  w).  The  2-D  inverse  Fourier  transform  of 
T{p,q,Zobj,TCs,oj)  gives  an  accurate  spatial  map  of  the  absorption  variations  6 ■,  y ,  Zobj)- 

The  reconstructed  images  from  experimental  data  are  shown  in  Fig. (2).  The  two  objects 
are  well  resolved  (with  a  peak-to-trough  of  >2).  The  reconstructed  x-y  positions  of  these  two 
slices  are  about  (-1.80,  -0.25)  cm  and  (1.85,  0.25)  cm,  close  to  their  true  x-y  positions  (-1.6, 
-0.3)  cm  and  (1.6,  0.3)  cm.  Inaccuracies  in  the  position  measurements  might  account  for  the 
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discrepancy.  Images  shown  in  Fig.  (2  (b-c))  are  unprocessed.  The  reconstructed  absorption 
coefficients  are  well  above  the  background  noise  level  and  close  to  the  true  values  (e.g., 
IJ'IT  !==^(0.22±0.03)  cm"^  )  and  nlf  rs(0.13±0.03)  cm“^.  The  noise  in  the  reconstruction  is 
mainly  from  the  finite  step  size  and  scan  region,  the  positional  error  and  the  electronics. 
The  refractive  index  mismatch  between  the  object  and  background  medium  also  contributes 
to  the  inaccuracy  in  the  reconstructed  optical  properties.  The  FFT  calculation  takes  less 
than  200  msec  CPU  time  on  Sun  SparclO  workstation  and  the  presence  of  multiple  objects 
does  not  increase  the  computation  complexity. 

The  feasibility  of  FFT  algorithm  for  imaging  scattering  object  has  also  been  demon¬ 
strated.  The  geometry  is  the  same  as  in  Fig.l  (a).  In  this  case  we  used  spherical  objects 
as  heterogeneities  instead  of  thin  slices  to  test  how  accurate  our  algorithm  is  for  imaging 
extended  objects.  Sphere  1  of  radius  0.75  cm  with  ^  16  cm"^  was  placed  at  (-1.6,  -0.3, 
3.0)  cm,  sphere  2  of  the  same  radius  with  «  26  cm"^  was  placed  at  (1.6,  0.3,  3.0)  cm. 
Both  spheres  have  the  same  absorption  coefficient  as  the  background,  e.g..  0.02  cm“h  We 
see  from  Fig. 3  that  positional  information  about  these  two  scattering  objects  is  recovered. 
In  this  case,  we  use  a  slice  through  the  sphere  center  at  z=Zobj  and  obtain  2-D  scattering 
contrast  image.  We  do  not  expect  to  reconstruct  the  scattering  coefficients  accurately  since 
the  objects  are  extended.  Notice  however  that  we  were  still  able  to  obtain  the  correct  po¬ 
sitions  and  we  see  that  the  two  objects  are  well  resolved  (Fig. 3).  Furthermore,  the  ratio  of 
the  reconstructed  scattering  coefficients  is  close  to  the  true  ratio,  .e.g.,  g's2  ~  1-35 

while  the  true  ratio  is  about  1.62. 

Depth  information  is  required  to  obtain  full  3-D  images  with  this  diffraction  tomography 
technique.  One  simple  method  is  to  use  a  secondary  localization  scheme  to  deduce  object 
depth.  Alternatively  two  projection  images  of  the  sample  along  orthogonal  directions  provide 
sufficient  information  for  3-D  reconstruction. 

To  conclude,  we  have  successfully  applied  near  field  diffraction  tomography,  diffuse  pho¬ 
ton  density  waves,  and  FFT’s  to  obtain  projection  images  of  hidden  objects  in  highly  scatter¬ 
ing  tissue  phantom.  It  may  be  possible  to  obtain  clinical  projection  images  in  real  time  with 
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this  fast  FFT  approach.  The  geometry  used  so  far  has  been  infinite.  In  practice,  boundary 
effects  present  important  problems.  On  the  one  hand,  matching  materials  might  be  used 
to  reduce  the  boundary  effects.  On  the  other  hand,  the  introduction  of  a  surface  integral 
term  in  Eq.(l)  or  better  Green  functions  (which  vanish  on  the  extrapolated  boundary)  may 
be  used  to  incorporate  boundary  effects.  The  technique  presented  provides  a  basis  for  more 
complicated  and  realistic  reconstruction  methods  to  address  these  issues. 

It  is  a  pleasure  to  acknowledge  useful  discussions  with  Maureen  O’Leary.  A.G.Y  acknowl¬ 
edges  support  from  the  NSF  under  grant  No.  DMR93-06814.  B.C.  acknowledges  support  in 
part  from  NIH  under  grant  Nos.  CA  50766  and  CA  60182. 
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FIG.  1.  (a)  and  (b)  illustrate  the  experimental  geometry  and  setup,  respectively.  In  (a):  the 
source  is  fixed  at  the  origin  and  the  detector  scans  in  a  planar  geometry  at  z=Zf^~5-0  cm  over  9.3 
X  9.3  cm^  region.  In  (b):  APD- Avalanche  Photo  Diode;  A-AmpMfier;  SSB  I/Q-Single  Side  Band 
In-phase/ Quadrature-phase  Demodulator:  LO-Local  Port;  RF-Radio  Frequency  Port;  ADC- Analog 
to  Digital  Converter;  GPIB-General  Purpose  Interface  Bus. 
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FIG.  2.  (a)  shows  exact  positions  of  the  two  slices,  (b)  and  (c)  illustrate  the  surface  and  planar 
images  of  two  absorbing  slices.  The  absorption  coefficients  and  the  positions  of  the  two  slices  are 
recovered  as  shown  in  (b)  and  (c). 
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FIG.  3.  (a)  shows  exact  positions  of  the  two  spheres,  (b)  and  (c)  illustrate  the  surface  and 
planar  images  of  two  scattering  spheres.  The  FFT  algorithm  is  sensitive  to  the  scattering  properties 
and  the  reconstructed  ratio  of  scattering  coefficients  is  close  to  the  true  ratio. 


